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Abstract 



In this paper the author discusses the article by the HERMES Collaboration, 
Phys. Lett. B 666, 446 where several important results concerning strange 
quark properties in the nucleon were presented. The strange sea distribution 
(]^ \ was found to be very different from the non-strange in the Bjorken x scaling 

• variable. In addition, the magnitude of these two distributions at low x is 

similar, contrary to most of the available Parton Distribution Function sets. 
The strange quark helicity distribution was also extracted combining inclusive 
and semi-inclusive kaon asymmetries. The result of the first moment of the 
strange quarks helicity distribution in the measured range is slightly positive, 
while a rather significant negative value is expected from the world polarised 
inclusive Deep Inelastic Scattering measurements. 

The author shows that a certain combination of fragmentation functions 
extracted from the preliminary HERMES kaon multiplicities presents a very 
strong dependence. Such a strong dependence is not expected from the 
DGLAP evolution of the fragmentation functions. If a similar dependence was 
present in the data used in the analyses of the discussed article, their conclusion 
might have to be altered. 
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1. Introduction 

The information on the strange quark properties in the nucleon is rather 
scarce. The strange quark density is poorly known, the main information comes 
from the neutrino experiments NuTeV 2] and CCFR [3]. The analyses 043 
suggest that the strange sea is suppressed with respect to the non-strange sea 
by a factor of about 2. This fact is not observed in the analysis of 8]. 

The HERMES results brought a very important contribution to the 
understanding of the strange quark properties in the nucleon. 
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First of all it was shown that the strange quark parton distribution function 
(PDF) has very different shape in the Bjorken x scaling variable as compared to 
the distributions of u and d. It was also verified that for x « 0.04 the densities 
of the strange and non-strange sea are comparable. This result was recently 
confirmed by the ATLAS analysis foll, but the final uncertainty of the ATLAS 
result is still large. 

Even more interesting results were obtained by HERMES in the polarised 
strange sector. From deep inelastic scattering measurements of the spin de- 
pendent structure function gi it is known that the first moment of the strange 
quark distribution is negative /\S = As(a;) + As{x)dx = -0.09 ± 0.01 ± 0.01, 
(see e.g. [13, [HI)- However, HERMES analysis based on the combination 



of the inclusive asymmetry and the semi-inclusive kaon asymmetry ^ 
concluded that in the measured range of x, the corresponding AS is consis- 
tent with zero and negative values are not preferred; Jq q2 ^^i^) + As{x)dx = 
0.037 ± 0.019 ± 0.027. These results are often called strange quark polarisation 
puzzle. 

The HERMES observation lead to a sign changing solution of the polarised 
strange quarks distribution function in a PDF set by D. de Florian, R. Sassot, 
M. Stratmann and W. Vogelsang 12|. In the high x region the strange quark 
polarisation is positive, changing sign for low values of x. As a result the semi- 
inclusive and the inclusive analyses of AS do no longer contradict each other. 

The other possible solution of the puzzle was proposed by E. Leader, A. V. 
Sidorov and D. B. Stamenov They argued that the puzzle only appears if 
one uses the DSS set of fragmentation functions (FFs) [3], as done in [l[. Using 
another set, namely the HKNS [ll| one, semi-inclusive data prefer negative 
strange polarisation in the whole x range. Based on published materials in [l| 
it can be shown that the proposition of [31 has a big weakness. Namely the 
two discussed sets of FFs differ indeed a lot in the kaon sector. The strange 
quark fragmentation function into kaons is more than a factor two higher in DSS 
than in HKNS. As a result if the HKNS FF had been used for the extraction 
of the unpolarised strange quark PDF from the HERMES data, the resulting 
distribution would have been more than a factor two higher than extracted in 

and in disagreement with all experiments. 

However, the results of [l|, published in 2008 results, are based on the anal- 
ysis of the sum of kaon (K^) multiplicities. The author analysed non-published 
but publicly available, preliminary HERMES results [16] presented in 2011, as 
well as results shown in PhD thesis of [l7j based on HERMES data, from 2005. 
In these two the sources charged separated kaon multiplicities are available. It 
will be shown that a certain combination of FFs extracted from the difference of 
multiplicities has a very strong (x) dependence not expected by DGLAP. 
If such a dependence was present in the data used in [l|, the explanation of E. 
Leader et al. could not be excluded, while conclusion of [l[ might had to be 
changed. 
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2. HERMES analysis of unpolarised strange sea 

The HERMES results [l[ on the strange quark PDF S(x) are based on the 
analysis of sum ofK^ multiplicities dN'^{x)ldN^'^{x) in semi-inclusive deep 
inelastic scattering (SIDIS) of electrons on a deuteron target. As stated there: 



S{x) / D^{z)dz ~ Q{x) 



d^N"^ (x) 



-J D^{z)dz 



(1) 



where S{x) — s{x) + s(x), z is the ratio of energies of the hadron and the 

K ( 



virtual photon in the target rest frame, Dq{z) = 4L)^g(z) + Djf Jz) and 



D^{z)^2D^^-M- . n 

As presented in Fig. 1 of the kaon multiplicity is flat for high x, i.e. in 
the region where there is no strange quarks and it rises by about 50% for lower 
values of x. Without strange quarks the distribution should be almost flat (or 
even decreasing for lower a;) therefore the large excess of kaon multiplicity in the 
low X region is interpreted as a strong signature of the strange quarks presence. 
Their contribution was quantified using Eq. ([T]). The value of J Dg {z)dz was 
taken from DSS and / Dq {z)dz was extracted directly from the HERMES data 
at high X. In the first order, the HERMES analysis neglected the possible 
negative four momentum transfer dependence of the FF. However, in the 
later stage of the analysis the dependence as in DSS fragmentation function 
is taken into account. The expected dependence is weak, of the order of 
10-20%, in the measured range of ^ (i _ iq) [GeV/cf. 



3. Tentative analysis of the HERMES publicly available data 

The author extracted the and K'^ multiplicities on deuteron target 
from slide 15 of [l^. These multiplicities cover the following kinematic ranges: 
X e (0.03 - 0.5) and G (1 - 10) (GeV/c)^, and they were integrated in the 
z range from 0.2 to 0.8. To be able to perform some quantitative analyses, 
the author of this paper assumes that in the strange sector there are three 
fragmentation functions: Dstr = = Df ; Djav = = D^^ and D„„/ 

for the remaining combinations. Thus: Dq = 4:D!^_^^ + -D^^ ~ 4:Dfav + GDunf 
and Dg = 2_D^j — 2Dstr + 2I?„„/. Here for simplicity it was assumed that 
Di = J°^^ Di{z)dz. 

It turned out that a very important distribution to understand better the 
strange sector is the difference of and multiplicities, dNf^^^ (x) / dN-^^^{x) . 
First of all the multiplicity difference leads to a non-singlet distribution, whose 
evolution decouples from gluons, and is thus simpler. Secondly a lot of sys- 
tematic uncertainties cancel in the multiplicity difference. Thirdly the disputed 
strange fragmentation function does not contribute to the difference. Finally 
using Eq. ([2]) one has a direct link to a certain combination of FF, namely 

D fav Dy^jiJ . 

dNf,ff{x) _ A{u,{x) + d,{x)) 



^^D/S(2,) ^[u{x) + u{x) + d{x) + d{x)) + As{x) 



{Dfav{x)-Dunf{x)) (2) 
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The presented equation is for the deuteron target. The unpolarised quark dis- 
tributions of various flavors are denoted by u{x), u{x), d{x), d{x), s{x). The de- 
nominator of the equation is closely related with the measured I7 cross-section, 
while Uy{x) and d^ix) are the valence distributions of the u and d quarks re- 
spectively. Therefore both these factors are known experimentally with a good 
precision. There is a strong correlation between x and in the HERMES 
experiment, thus for simplicity it was assumed here that FFs depend upon x, 
instead of • 

The author evaluated Dfav{x) — Dunf{x) using Eq. ^ and the CTEQ6L 



PDF set [18| . The systematic errors of and K~ multiplicities are strongly 
correlated; but the correlation factors were not known to the author. He decided 
to scale the errors as given in Q so that the fit of {Dfav — Dunf){logio{x)) 
by a straight line gives the value equal to the number of degrees of freedom. 
The scaling factor was found to be 2.5, thus quite large. 

The resulting Dfav{x) — (x) are presented in Fig.[T] The points are the 
Dfav{x) — Dun fix) values integrated in the measured range of z G (0.2 — 0.8), 
extracted by the author. The solid line is the straight line fit to these data. 
Finally the dashed-line is the DSS parametrisation, taking into account the 
average values in the different x intervals. The values of the multiplicities 
and PDF parameters used by the author are summarised in Tables [T] and [5J 

Contrary to the expectation of 1J|, a very strong x or dependence of 
the evaluated FFs is observed. The author also verified that a similar, but a 
factor 2.5 weaker, dependence is also present in the results from the PhD thesis 
of However, for those results the multiplicities were evolved to the same 

value, and in addition a region in z smaller than in [l[ was available. The 
author of [l3| mentions that there is possibly a problem in the procedure of the 
multiplicity extraction, since the extracted K~ multiplicities were sometimes 
negative. 



"°^0.12 

0.08 
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0.04 
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Figure 1: The distribution of f^'^ Dj^^ — D^„fdz as a function of x, obtained by the author 
from HERMES prel. data [3]. The continuous line is a straight line fit, while the dashed line 
corresponds to the DSS predictions related with dependence of fragmentation functions. 



HERMES used the high x region to constrain Dq , keeping this value for 
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low X or slightly changing it to take into account weak evolution from DSS. 
However, this assumption is not supported by their preliminary data where 
evolution of FFs can be studied. Since Dq = 4(£)/ati — -D„„/) + 10-D„„/, a strong 
X or Q"^ dependence of D fav ~ influences • The possible impact of 

the observed x or dependence on the results obtained in [ll] is discussed in 
two scenarios. 

In the first scenario Dunf is large at high x and rapidly decreases for lower x 
values. A typical in DSS is of the order of 0.012-0.015; However, the ob- 

served change of D fav (x) — D^nf (x) is a factor about 5 larger. Thus introducing 
a strong x or dependence of Dunf results in negative values in most of the 
phase space. At the same time the strange quark density is further increased. 

In the second scenario Dfav, *.e. and , increases for low x. In 

this case the increase of the sum of K'^, K~ multiplicities for low x observed 
in HERMES corresponds mainly to kaons originating from u, u instead of s, s. 
As a result the strange quark PDF is lower than anticipated in [l[. It follows 
naturally that the non-negative AS* obtained in the HERMES analysis could be 
explained by the fact that more than expected, positively polarised, u, u quarks 
contribute to a kaon sample from which the spin dependent asymmetry ^fT^ is 
extracted. As a consequence the obtained AS* is biased towards positive values. 

If the discussed data were used for the extraction of FFs, in the discussed 
scenario the resulting Dgtr would be overestimated. As in previous cases part 
of the kaons originating from u and u would be wrongly assigned as having s 
or s origin. This possible overestimation of Dgtr has further infiuence on the 
A5' extraction, resulting in too optimistic statistical uncertainties and possibly 
further bias towards a positive value. An example of the correlation of the 
uncertainty and the value of AS* for different ratios of Dstr/ Dfav can be seen 
in COMPASS analysis 19]. 

The author tried to verify a possible impact of the observed Dfav{x) — 
Dunf{x) dependence on the Dstr- In order to do so a simple LO extraction of 
FFs from HERMES data was performed using as input the extracted integrated 
K+ and R- muhiphcities from [11] and the CTEQ6L PDF set. 

Without any Q"^ dependence of FFs the fit ^ jndj was 75.4/15. The fol- 
lowing z-integrated values of FFs were obtained: D^tr = 0.45 ± 0.09, Dfav = 
0.100 ±0.003 and = 0.017 ±0.002. For comparison, the results of LO DSS 

analysis of the world data at Q^— 2.5 (GeV/c)^ are presented: Dstr = 0.62, 
Dfav = 0.091 and D„„/ = 0.0120 . In the next step, an additional free param- 
eter was added, modifying Dfav{x) = Dfav{^) + par ■ logiQ{x). The functional 
form of the x dependence is selected taking in to account on the empirical ob- 



^ Observe, that this conclusion is valid even in case one assumes only charge conjunction 
asymmetry of FFs, thus 6 independent FFs in the strange quark sector. 

■^As a matter of fact some of the multiplicities obtained from HERMES data were used in 
the DSS FF fit, namely Information from [T?! together with private communication with the 
author of that PhD thesis. However, detailed analysis of the impact of HERMES data and 
their unforeseen Df^^ix) — D^^f{x) dependence on the full FF fit is far too complex to be 
done by the author in this paper. 
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servation made in Fig. [TJ The resulting of the fit improved by a factor 4 to 
18.1/14. The slope was found to be about -0.061 as can be deduced from Fig.[T]. 
However, at the same time the strange FF changed a lot, Dstr = —0.44 ± 0.14, 
i.e. the central value became unphysical. Adding an additional dependence 
for Dunf as expected from the DSS brings the central value of Dgtr to -0.3 or 
-0.2 depending on if the evolution of FF from DSS is applied as an additive 
or a multiplicative correction. Thus the allowed region is concentrated only for 



low values of Dstr «-e. as obtained in [15| and discussed in [13 



For completeness the author also verified that if the D fav is kept constant 
while Dunf lias a x dependence of the form Dunf{x) = D„„/(l) +par ■ logio{x), 
the resulting x^/ndf is about 37.5/14, but Dunf is negative for x < 0.15 and 
in addition Dstr/Dfav ~ 25 thus not very reliable. A simultaneous fit of the x 
dependence of both Dfav and D„„/ did not improve the /ndf value signif- 
icantly: 17.6/13 to be compared with 18.1/14 when only Dfav{x) dependence 
was included. In addition the correlation between the fit parameters exceeds 
99% in this case. 

However, the most important observation is that indeed there is a very 
strong correlation between Dfav{x) — Dunf(x) and Dstr- Without a proper 
understanding of the x dependence of D fav — Dunf it is very hard to obtain from 
SIDIS data any non-disputable results in the strange quark sector. Observe that 
a systematic effect giving rise to a 10% slope of Dfav{x) — Dunf{x) results in a 
Dstr/Dfav bias of about 1 . The author doesn't speculate why such a strong x or 

dependencies are observed in the preliminary HERMES data sets. However, 
it is worth to mention that for the preliminary COMPASS multiplicities [20l |. 
the author extracted the Dfav{x) — Dunf{x), which within uncertainties is flat 
in the whole measured x region x £ (0.004 — 0.7) . The range is between 
1-30 (GeV/c)2 for these data. 



4. Summary 

Based on publicly available HERMES data the author showed that the anal- 
ysis of K'^ and K~ multiplicities results in a not foreseen x dependence of the 
Dfav{x) — Dunf(x) fragmentation functions. This dependence can explain, in 
a certain scenario, the different than expected strange quark PDF obtained by 
HERMES as well as it will bias the strange quark helicity distribution towards 
a positive value. Therefore the key point is the verification by the HERMES 
collaboration that D fav — Dunf does not show any unforeseen x dependence for 
the data used in [l[. The answer to this question is very important in order to 
understand the strange quark sector in the fixed target experiments domain. 
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Appendix 



In order to make the verification of the presented results possible, in Tables 
[T] and [2] the raw information used in this analysis is given. 



Table 1: Mean values of x and in (GeV/c)^ for the used data as well as values of the 
a;-PDF for quark flavours from CTEQ6L. In the last column 4(mj, + d^) / {h{u + u + d + d) + is) 
is given. 



X 




u 


u 


d 


d 


s 


4(ti„+d„) 
b{u+u+d+d)+As 


0.033 


1.1 


0.408 


0.165 


0.346 


0.178 


0.063 


0.286 


0.048 


1.4 


0.450 


0.153 


0.359 


0.178 


0.062 


0.322 


0.065 


1.6 


0.488 


0.138 


0.367 


0.176 


0.057 


0.357 


0.087 


1.7 


0.527 


0.119 


0.372 


0.169 


0.050 


0.399 


0.120 


2.1 


0.579 


0.093 


0.375 


0.150 


0.040 


0.462 


0.170 


3.1 


0.624 


0.064 


0.358 


0.108 


0.027 


0.551 


0.240 


4.9 


0.618 


0.039 


0.302 


0.053 


0.015 


0.646 


0.340 


7.4 


0.513 


0.019 


0.202 


0.014 


0.006 


0.727 


0.450 


10.1 


0.344 


0.006 


0.104 


0.002 


0.002 


0.767 



Table 2: — K~ multiplicity differences extracted from protect [T^ and their errors. Differ- 
ent intervals of x and z are shown. The final 2-integrated multiplicity difference is presented 
in the last column. 



X 




z e (0.2 - 0.3) 


z e (0.3 - 0.4) 


z e (0.4-0.6) 


z e (0.6-0.8) 


dNOis(x) 


0.033 


1.1 


0.084 ±0.026 


0.080 ±0.023 


0.056 ±0.009 


0.031 ±0.008 


0.0338 ±0.0042 


0.048 


1.4 


0.084 ±0.012 


0.095 ±0.016 


0.054 ±0.006 


0.032 ±0.006 


0.0352 ±0.0026 


0.065 


1.6 


0.090 ±0.023 


0.084 ±0.007 


0.061 ±0.007 


0.035 ±0.004 


0.0367 ±0.0029 


0.087 


1.7 


0.088 ±0.035 


0.081 ±0.017 


0.071 ±0.005 


0.039 ±0.009 


0.0391 ±0.0044 


0.120 


2.1 


0.097 ±0.041 


0.082 ±0.014 


0.067 ±0.006 


0.040 ±0.003 


0.0395 ±0.0046 


0.170 


3.1 


0.091 ±0.052 


0.091 ±0.016 


0.074 ±0.005 


0.037 ±0.006 


0.0404 ±0.0057 


0.240 


4.9 


0.116 ±0.027 


0.096 ±0.011 


0.078 ±0.006 


0.041 ±0.004 


0.0451 ±0.0033 


0.340 


7.4 


0.107 ±0.045 


0.108 ±0.016 


0.087 ±0.008 


0.040 ±0.003 


0.0469 ±0.0051 


0.450 


10.1 


0.117 ±0.088 


0.028 ±0.048 


0.071 ±0.008 


0.027 ±0.004 


0.0341 ±0.0102 
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